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We propose to use diatomic molecular ions to search for strongly enhanced effects of variation of 
fundamental constants. The relative enhancement occurs in transitions between nearly degenerate 
levels of different nature. Since the trapping techniques for molecular ions have already been devel¬ 
oped, the molecules HBr + , HI + , Br]", I]", IBr + , IC1 + , and IF + are very promising candidates for 
such future studies. 

PACS numbers: 06.20.Jr, 33.20.Vq, 33.15.Mt, 33.15.Pw, 33.57.+C 


Theories unifying gravity with the other interactions 
suggest the possibility of spatial and temporal variation 
of dimensionless fundamental physical constants, such as 
the fine structure constant, a = e 2 /he and the proton 
to electron mass ratio, /j, = m p /m e mm, the latter 
being connected to the strong coupling constant Oil!. 
Two major directions for the search of a and t 1 varia¬ 
tion are observational studies, such as the analyses of 
high-resolution quasar absorption spectra (Refs. [5H?] 
and references therein), and laboratory research, includ¬ 
ing atomic and molecular clocks 011- 

Molecules are very promising probes for variation of 
fundamental constants (VFC) as molecular spectra can 
be very sensitive to both a and /i Gann, making it pos¬ 
sible to look for variation of both constants in a single 
experiment. So far, the use of molecules has mostly fo¬ 
cused on astrophysical observations of molecular spectra 
at high red-shifts. Molecular probes used in astrophysical 
research include the Lyman and Werner absorption lines 
in H 2 and HD molecules nam, the 18 cm 1 transition 
of OH pT]; the inversion spectrum of ammonia p~5l - ll7] . 
A—doubling transitions of CH [18], and the internal ro¬ 
tational transitions of interstellar methanol m- 

The most stringent laboratory molecular limit to date 
on the variation of /i is /i//r = (—3.8 ± 5.6) x 10 -14 yr _1 
[20] , set in a two-year-long experiment over the course of 
which vibrational transitions in the SF 6 molecule were 
compared with the Cs standard. So far, this is the only 
laboratory experiment using molecular systems to study 
VFC. However, the development of methods for creation 
and trapping of ultracold molecules [21] [22! opens ex¬ 
citing new prospects to test for VFC by means of high- 
resolution molecular spectroscopy. In such experiments it 
is important to use molecules that experience enhanced 
sensitivity to VFC. However, locating such favourable 
systems is not a trivial task. 

Nearly degenerate levels with different sensitivity to 
VFC may provide huge enhancements of a relative vari¬ 


ation, since 5lo/ui tends to infinity when the distance be¬ 
tween the levels ui is close to zero. This is the case with 
dysprosium 0123 . Relatively low accuracy of the transi¬ 
tion frequency measurements in Dy provides record-high 
accuracy in VFC measurements [23]. However, the prob¬ 
lem with Dy is that one of the levels has a short lifetime, 
thus limiting the accuracy. 

Molecules can have both degenerate levels and very 
narrow transition widths. Diatomic molecules that have 
a near cancellation between hyperfine structure and rota¬ 
tional intervals [25] or between fine structure and vibra¬ 
tional intervals [26] can benefit from such an enhanced 
sensitivity. A number of such molecules have already 
been identified, e.g. CS 2 [27] [25], CaH, MgH, CaH+ 
J2H3U, Cl]", IrC, HfF + , SiBr, LaS, LuO, and others 
[25l [251 [52] . Several polyatomic molecules were suggested 
for laboratory studies of VFC as well, see e.g. [33] EH 
and references in M- Here we note that practical 
advantage can be gained by using molecular ions rather 
than neutral molecules due to the greater ease of trap¬ 
ping and cooling [351 136] . Unfortunately, experimental 
spectroscopic constants for charged diatomics are scarce 
compared to their neutral counterparts, and the search 
for VFC using such molecules is somewhat limited. 

In this work we investigate a group of diatomic 
molecules that satisfy all of the above mentioned crite¬ 
ria, namely the singly charged dihalogens and hydrogen 
halides. We are interested in the forbidden rovibrational 
transitions between the nearly degenerate sub-levels of 
the X 2 U ground states of these molecules. Using avail¬ 
able experimental spectroscopic data, we identify such 
promising transitions, and examine the dependence of 
the transition energy on both a and /r. Strong enhance¬ 
ment of the sensitivity of these transitions to variations 
in a and in /1 makes them excellent candidates for future 
experiments to detect VFC. 

The systems presented here were selected for an in- 
depth study due to the availability of their accurate 


2 


experimental spectroscopic parameters, which are listed 
(along with their definitions) in the Supplementary Ma¬ 
terial at [URL will be inserted by publisher]. These 
systems include four isotopologues of HBr + (H 79 Br + , 
B sl Br+, D 79 Br+, and D 81 Br+) gH EE], two isotopo¬ 
logues of BI + (BI + and DI + ) EH], three isotopologues 
of Br+ ( 79 Br+, 81 Br+, and 79 Br 81 Br+) gOJ, 1+ gUlll], 
two isotopologues of IBr + (I 79 Br + and I 81 Br + ) [43], two 
isotopologues of IC1+ (I 35 C1 + and I 37 C1 + ) [M], and IF+ 
gS]. A number of other charged dihalogens and hydro¬ 
gen halides were considered, but were found unsuitable 
for the present study. These systems are discussed in 
more detail in the Supplementary Material. 

Using available spectroscopic constants, we reproduce 
the molecular potential energy curves and the spectro¬ 
scopic levels by following the Rydberg-Klein-Rees pro¬ 
cedure gHHH]. This allows us to locate the promising 
low-energy rovibrational transitions between the two sub¬ 
states of the X 2 H ground state of the molecules of inter¬ 
est. Such low energy transitions will exhibit a strong 
enhancement of sensitivity to VFC. Once the favourable 
transitions are located, we examine the dependence of 
their energy on the a and y, constants. Following the 
line of thought from earlier work gB] EJl we derive sim¬ 
ple analytical expressions for the transition energies and 
their fractional variation. Rere we do not deal with the 
lowest vibrational and rotational levels only, but also in¬ 
clude the effect of anharmonicity and rotational spacings, 
which was neglected in earlier publications. Taking only 
the harmonic frequency w e , the first anharmonicity u} e x e , 
the rotational constant B e , and the spin-orbit coupling 
constant A e and its first vibrational dependence into 
account (thus neglecting higher order terms), the transi¬ 
tion energy to between the two states can be expressed 
as 


u> = E' — E" = w e (z/ - v") 

— uj e x e (y' — i'")(y" + v' + 1) 

+ B e {J'- + J'+ 1) (1) 

-A e -^AW(v" + v' + 1 ), 

where we denote the initial state by double prime and the 
final state by single prime; in this case these correspond 
to the X 2 Ils and the A 2 Ri states, respectively. 

The spin-orbit coupling constant A v and its con¬ 
stituents A e , A^\ A( 2 \ etc. scale as ~ a 2 . The rest of 
the spectroscopic constants are only weakly dependent on 
a and this dependence can be neglected [32]. Following 
from the dependence of various spectroscopic constants 
on the reduced mass [5B], we can easily deduce their scal¬ 
ing with the proton to electron mass ratio /r: u> e and A^ 
scale as ~ M _I , while u> e x e and B e scale as ~ /r -1 . 

Using these analytical dependencies along with Eq. |T] 
the fractional variation of the transition energy to can be 



FIG. 1. (color online) Potential energy curves of the two 
substates of the X J n state in H' 9 Br + . Arrows point to the 
quasi-degenerate vibrational levels. Origin of the vertical axis 
corresponds to the bottom of the unsplit potential well. 


calculated as 


Sui 


5fi 


5a 


— =K^+K a — 


( 2 ) 


with the relative enhancement factors K M and K a defined 
as 


K M =[^oj e (iy" - i/) - uj e x e (y" - v')(v" + v' + 1) 

+ B e (J" - + J'+ 1) (3) 

+ iA (1) (j/' + v' + l)]w -1 = A^w -1 , 


K a = [~2A e - AW(v" + i/ + l)]^- 1 = KaU- 1 , (4) 

where the symbols K a and K M represent the absolute 
enhancement factors. 

We have examined the effect of including the next or¬ 
der parameters in the Dunham series, i.e. uj e y e , a e , and 
A( 2 ), on the calculated K^ and K a , and found it to be 
negligible for the majority of the cases considered here. 

We illustrate the scheme outlined above by using 
R 79 Br+ as an example. Figure [l] presents the recon¬ 
structed potential energy curves of the two sub-states 
of the X 2 R ground state of R 79 Br + . We find that the 
Af 2 U|, v = 14 level is very close in energy to the X 2 H i, 
v = 12 level, differing by 19 cm -1 . We then take a 
closer look at the rotational states of the favourable vi¬ 
brational transition, in order to select even closer quasi¬ 
degeneracies of the energy levels (Figure [2] left). Vibra¬ 
tional levels V 2 R3, v = 15 and A' 2 Ri, v = 13, although 
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FIG. 2. (color online) Rotational structure of the selected 
quasi-degenerate vibrational levels of H 7!) Br + with transitions 
represented by dashed arrows. Allowed e e and / / 

transitions are denoted with e and /, respectively, in the 
brackets with transition energies (cm -1 ). 


futher apart (123 cm -1 ), also exhibit quasi-degeneracies 
in the rotational structure (Figure |2j right). 

We focus only on the dominant single-photon dipole 
transitions obeying the A J = 0, ±1 selection rule. Since 
the transitions conserving the total angular momentum 
(A J = 0) can be blended together for all values of J 
[52] . thus negatively affecting the accuracy, we are inter¬ 
ested in the A J = ±1 transitions only. Furthermore, we 
also include the effect of A-doubling and we only consider 
transitions following the rotational symmetry selection 
rule e<«/. The A-doubling of the 2 Il3 rotational levels 
is negligible and not even visible in the scope of Figure 
[2j We can identify a number of rovibrational transitions 
corresponding to a unit change in the total angular mo¬ 
mentum J while having very low transition energies be¬ 
tween the reconstructed rovibrational levels. These are 
listed in Table [TJ along with the transition energies and 
the corresponding enhancement factors as defined in Eqs. 
[3] and [TJ The same procedure is performed for the rest 
of the systems under study. These results are also listed 
in Table [TJ and the corresponding plots are presented in 
the Supplementary Material. Note that there are more 
favorable transitions identified in the rotational plots, as 
these usually occur in series for several neighbouring ro¬ 
tational levels, however, we only present the ones with 
lowest ui (and therefore highest K M and K a ) in Table [I] 

The obtained values of K a look very promising and 
span orders of magnitude ranging from 10 3 to 10 6 . Sev¬ 
eral transitions in ij, I 81 Br + , and I 35 C1 + are particularly 
sensitive to variation in a, with K a values larger than 


TABLE I. Transition energies lo (cm 1 ) and corresponding 
enhancement factors, absolute, K M and K a (cm -1 ), and rel¬ 
ative, Kfj, and K a , for several identified transitions of cations 
under investigation. Symbols e and / denote the transitions 
of type ef> e and /«-»/, respectively. 


cation 

v" 

t/ 

J" 

J' 

S LJ 



K a 

K a 

H' y Br+ 

14 

12 

1.5 

2.5 

f 0.74 

-80.8 

-110 

5250 

7120 




4.5 

5.5 

/ 1-26 

129 

-103 


4180 


15 

13 

11.5 

12.5 

/ 0.14 

-426 

-3030 

5250 

37200 

H 81 Br+ 

14 

12 

1.5 

2.5 

/ 0.82 

-80.4 

-97.9 

5250 

6390 




5.5 

6.5 

/ -0.27 

-145 

543 


-19700 


15 

13 

11.5 

12.5 

e 1.44 

-426 

295 

5250 

-3630 




6.5 

5.5 

/ 0.52 

119 

-228 


10100 

D 79 Br+ 

9 

7 

7.5 

8.5 

/ 2.20 

881 

400 

5270 

2390 


10 

8 

1.5 

0.5 

e 1.99 

871 

437 

5270 

2650 

D 81 Br+ 

9 

7 

7.5 

8.5 

/ 2.49 

881 

354 

5270 

2120 


10 

8 

1.5 

0.5 

e 2.23 

871 

391 

5270 

2370 

HI+ 

6 

3 

3.5 

4.5 

e -3.05 

2060 

-675 

10700 

-3490 




21.5 

20.5 

/ -0.41 

2390 

-5780 


-25800 




26.5 

25.5 

e 2.63 

2450 

934 


4060 


12 

8 

7.5 

8.5 

/ 1.50 

984 

656 

10600 

7070 




21.5 

20.5 

e -2.35 

1360 

-580 


-4510 


16 

11 

16.5 

15.5 

e 0.12 

190 

1630 

10600 

91000 

DI+ 

7 

3 

1.5 

2.5 

e 0.27 

2230 

8220 

10700 

39300 




21.5 

20.5 

/ 0.98 

2390 

2420 


10800 




27.5 

26.5 

e 2.42 

2430 

1000 


4420 


14 

9 

9.5 

10.5 

e -0.91 

1460 

-1600 

10600 

-11700 




22.5 

21.5 

/ -0.76 

1670 

-2190 


-14000 




25.5 

24.5 

e -3.92 

1690 

-430 


-2710 


19 

3 

9.5 

10.5 

/ 1.08 

686 

638 

10600 

9860 




19.5 

18.5 

/ -3.64 

878 

-241 


-2910 

4 

30 

6 

30.5 

31.5 

e -0.00 

2270 

-682000 

9990 

-3010000 




34.5 

35.5 

/ -0.06 

2260 

-40000 


-176000 


39 

14 

50.5 

51.5 

e 0.01 

2090 

230000 

9800 

1080000 




54.5 

55.5 

/ -0.13 

2090 

-16200 


-76200 


48 

22 

1.5 

2.5 

/ 0.06 

1900 

30500 

9610 

155000 




19.5 

18.5 

e -0.02 

1900 

-79100 


-401000 




28.5 

27.5 

/ -0.04 

1900 

-45900 


-232000 

79 Br2+ 

11 

3 

22.5 

23.5 

e 0.15 

1350 

9170 

5570 

37900 




24.5 

25.5 

/ - 0.10 

1350 

-13700 


-56600 


31 

22 

18.5 

19.5 

e -0.10 

1190 

-12300 

5370 

-55400 




21.5 

22.5 

/ 0.04 

1190 

33200 


150000 

79 ’ 81 Br+ 

9 

1 

20.5 

21.5 

e 0.09 

1360 

14500 

5600 

59700 




22.5 

23.5 

/ - 0.11 

1360 

-12700 


-52200 


30 

21 

3.5 

4.5 

e -0.09 

1160 

-12900 

5450 

-60400 




4.5 

5.5 

/ -0.07 

1160 

-15500 


-73000 

00 

to 

8 

0 

12.5 

13.5 

e 0.11 

1370 

12900 

5600 

52700 




14.5 

15.5 

/ 0.09 

1370 

14500 


59400 


28 

19 

36.5 

37.5 

e 0.05 

1190 

25700 

5450 

118000 




39.5 

40.5 

/ - 0.10 

1190 

-11600 


-53400 

I 79 Br+ 

25 

8 

27.5 

28.5 

e -0.05 

1720 

-32200 

9220 

-173000 




30.5 

31.5 

/ 0.20 

1720 

8440 


45300 


39 

22 

28.5 

29.5 

e 0.04 

1020 

25200 

9140 

225000 




30.5 

31.5 

/ 0.02 

1020 

63400 


566000 

I 81 Br+ 

24 

7 

31.5 

32.5 

e 0.06 

1760 

31700 

9320 

168000 




34.5 

35.5 

/ 0.23 

1760 

7630 


40500 


36 

19 

7.5 

6.5 

/ - 0.01 

1170 

-83400 

9320 

-664000 




12.5 

11.5 

/ - 0.01 

1170 

-88700 


-707000 

I 35 C1+ 

20 

8 

18.5 

19.5 

e 0.31 

1940 

6160 

9550 

30300 




20.5 

21.5 

/ - 0.01 

1940 

-166000 


-817000 

I 37 C1+ 

20 

8 

18.5 

19.5 

e -0.11 

1940 

-17700 

9550 

-87000 




21.5 

22.5 

/ 0.15 

1940 

12800 


62900 

IF+ 

14 

5 

23.5 

24.5 

e -0.06 

2760 

-45900 

11200 

-186000 




26.5 

27.5 

/ -0.12 

2760 

-23500 


-95200 


22 

12 

34.5 

35.5 

e 0.49 

2780 

5680 

11100 

22600 




37.5 

38.5 

/ 0.30 

2780 

9260 


36900 


28 

17 

34.5 

35.5 

e -0.66 

2820 

-4270 

11000 

-16600 




37.5 

38.5 

/ -0.41 

2820 

-6820 


-26500 
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600,000, due to the almost perfect degeneracy between 
these levels, differing by only about 0.01 cm -1 . 

The calculated /\ M are generally one or two orders of 
magnitude smaller than K a . However, many theories 
suggest that the variation in /r should be on the order 
of 30 times larger than that in a (see e.g. [51] and ref¬ 
erences therein). Thus, the effect of possible fi -variation 
on the spectra of these systems should not be neglected. 
Since the first term in Eq. [3] is dominant, we can eas¬ 
ily see, that the larger the difference between initial and 
final state vibrational quantum number is, the higher is 
the enhancement factor A' M . This is automatically satis¬ 
fied for heavy diatomics with larger spin-orbit splitting 
between the two X 2 H sub-states as is illustrated by the 
results in Table HJ 

The identified transitions in the HBr + and HI + iso- 
topologues, as well as in 79 Br 81 Br + , may be of particular 
interest for experimental investigation because of their 
low- J nature. Rotationally cooled molecules should pro¬ 
duce spectra with higher resolution [52l - i55] . which is of 
great importance for detection of the minute changes in 
transition energies resulting from VFC. 

An advantage can be gained by considering the varia¬ 
tion of the ratio of two transition frequencies within the 
same rovibronic spectrum, W 1 /W 2 , rather than the varia¬ 
tion of the frequencies themselves [32]. The variation of 
uj\/u >2 is given by 


5{ui/w2) 

(wi/w 2 ) 


Sfj, 


5a 




a 


(5) 


Considering the variation of a dimensionless ratio of 
two transition energies removes the dependence on the 
unit system and minimizes the systematic effects. In ad¬ 
dition, selecting two transitions where either A' 7 and K 2 , 
or and K 2 , or both pairs, have an opposite sign, max¬ 
imizes the sensitivity of the measurement. 
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